Effects of repeated alcohol abstinence on within-subject prefrontal cortical gene expression in rhesus macaques.

Male rhesus monkeys (n = 24) had a biopsy of prefrontal cortical area 46 prior to chronic ethanol self-administration (n = 17) or caloric control (n = 7). Fourteen months of daily self-administration (water vs. 4% alcohol, 22 h access/day termed "open-access") was followed by two cycles of prolonged abstinence (5 weeks) each followed by 3 months of open-access alcohol and a final abstinence followed by necropsy. At necropsy, a biopsy of Area 46, contralateral to the original biopsy, was obtained. Gene expression data (RNA-Seq) were collected comparing biopsy/necropsy samples. Monkeys were categorized by drinking status during the final post-abstinent drinking phase as light (LD), binge (BD), heavy (HD) and very heavy (VHD drinkers). Comparing pre-ethanol to post-abstinent biopsies, four animals that converted from HD to VHD status had significant ontology enrichments in downregulated genes (necropsy minus biopsy n = 286) that included immune response (FDR < 9 × 10-7) and plasma membrane changes (FDR < 1 × 10-7). Genes in the immune response category included IL16 and 18, CCR1, B2M, TLR3, 6 and 7, SP2 and CX3CR1. Upregulated genes (N = 388) were particularly enriched in genes associated with the negative regulation of MAP kinase activity (FDR < 3 × 10-5), including DUSP 1, 4, 5, 6 and 18, SPRY 2, 3, and 4, SPRED2, BMP4 and RGS2. Overall, these data illustrate the power of the NHP model and the within-subject design of genomic changes due to alcohol and suggest new targets for treating severe escalated drinking following repeated alcohol abstinence attempts.

Undiagnosed RASopathies in infertile men.

RASopathies are syndromes caused by congenital defects in the Ras/mitogen-activated protein kinase (MAPK) pathway genes, with a population prevalence of 1 in 1,000. Patients are typically identified in childhood based on diverse characteristic features, including cryptorchidism (CR) in >50% of affected men. As CR predisposes to spermatogenic failure (SPGF; total sperm count per ejaculate 0-39 million), we hypothesized that men seeking infertility management include cases with undiagnosed RASopathies. Likely pathogenic or pathogenic (LP/P) variants in 22 RASopathy-linked genes were screened in 521 idiopathic SPGF patients (including 155 CR cases) and 323 normozoospermic controls using exome sequencing. All 844 men were recruited to the ESTonian ANDrology (ESTAND) cohort and underwent identical andrological phenotyping. RASopathy-specific variant interpretation guidelines were used for pathogenicity assessment. LP/P variants were identified in PTPN11 (two), SOS1 (three), SOS2 (one), LZTR1 (one), SPRED1 (one), NF1 (one), and MAP2K1 (one). The findings affected six of 155 cases with CR and SPGF, three of 366 men with SPGF only, and one (of 323) normozoospermic subfertile man. The subgroup "CR and SPGF" had over 13-fold enrichment of findings compared to controls (3.9% vs. 0.3%; Fisher's exact test, p = 5.5 × 10-3). All ESTAND subjects with LP/P variants in the Ras/MAPK pathway genes presented congenital genitourinary anomalies, skeletal and joint conditions, and other RASopathy-linked health concerns. Rare forms of malignancies (schwannomatosis and pancreatic and testicular cancer) were reported on four occasions. The Genetics of Male Infertility Initiative (GEMINI) cohort (1,416 SPGF cases and 317 fertile men) was used to validate the outcome. LP/P variants in PTPN11 (three), LZTR1 (three), and MRAS (one) were identified in six SPGF cases (including 4/31 GEMINI cases with CR) and one normozoospermic man. Undiagnosed RASopathies were detected in total for 17 ESTAND and GEMINI subjects, 15 SPGF patients (10 with CR), and two fertile men. Affected RASopathy genes showed high expression in spermatogenic and testicular somatic cells. In conclusion, congenital defects in the Ras/MAPK pathway genes represent a new congenital etiology of syndromic male infertility. Undiagnosed RASopathies were especially enriched among patients with a history of cryptorchidism. Given the relationship between RASopathies and other conditions, infertile men found to have this molecular diagnosis should be evaluated for known RASopathy-linked health concerns, including specific rare malignancies.

The human infertility single-cell testis atlas (HISTA): an interactive molecular scRNA-Seq reference of the human testis.

BACKGROUND: Single-cell RNA-seq (scRNA-Seq) has been widely adopted to study gene expression of the human testis. Several datasets of scRNA-Seq from human testis have been generated from different groups processed with different informatics pipelines. An integrated atlas of scRNA-Seq expression constructed from multiple donors, developmental ages, and fertility states would be widely useful for the testis research community. OBJECTIVE: To describe the generation and use of the human infertility single-cell testis atlas (HISTA), an interactive web tool for understanding human spermatogenesis through scRNA-Seq analysis. METHODS: We obtained scRNA-Seq datasets derived from 12 donors, including healthy adult controls, juveniles, and several infertility cases, and reprocessed these data using methods to remove batch effects. Using Shiny, an open-source environment for data visualization, we created numerous interactive tools for exploring the data, some of which support simple statistical hypothesis testing. We used the resulting HISTA browser and its underlying data to demonstrate HISTA's value for testis researchers. RESULTS: A primary application of HISTA is to search by a single gene or a set of genes; thus, we present various analyses that quantify and visualize gene expression across the testis cells and pathology. HISTA also contains machine-learning-derived gene modules ("components") that capture the entire transcriptional landscape of the testis tissue. We show how the use of these components can simplify the highly complex data in HISTA and assist with the interpretation of genes with unknown functions. Finally, we demonstrate the diverse ways HISTA can be used for new data analysis, including hypothesis testing. DISCUSSION AND CONCLUSIONS: HISTA is a research environment that can help scientists organize and understand the high-dimensional transcriptional landscape of the human testis. HISTA has already contributed to published testis research and can be updated as needed with input from the research community or downloaded and modified for individual needs.

TAD evolutionary and functional characterization reveals diversity in mammalian TAD boundary properties and function.

Topological associating domains (TADs) are self-interacting genomic units crucial for shaping gene regulation patterns. Despite their importance, the extent of their evolutionary conservation and its functional implications remain largely unknown. In this study, we generate Hi-C and ChIP-seq data and compare TAD organization across four primate and four rodent species and characterize the genetic and epigenetic properties of TAD boundaries in correspondence to their evolutionary conservation. We find 14% of all human TAD boundaries to be shared among all eight species (ultraconserved), while 15% are human-specific. Ultraconserved TAD boundaries have stronger insulation strength, CTCF binding, and enrichment of older retrotransposons compared to species-specific boundaries. CRISPR-Cas9 knockouts of an ultraconserved boundary in a mouse model lead to tissue-specific gene expression changes and morphological phenotypes. Deletion of a human-specific boundary near the autism-related AUTS2 gene results in the upregulation of this gene in neurons. Overall, our study provides pertinent TAD boundary evolutionary conservation annotations and showcases the functional importance of TAD evolution.

Consensus label propagation with graph convolutional networks for single-cell RNA sequencing cell type annotation.

MOTIVATION: Single-cell RNA sequencing (scRNA-seq) data, annotated by cell type, is useful in a variety of downstream biological applications, such as profiling gene expression at the single-cell level. However, manually assigning these annotations with known marker genes is both time-consuming and subjective. RESULTS: We present a Graph Convolutional Network (GCN)-based approach to automate the annotation process. Our process builds upon existing labeling approaches, using state-of-the-art tools to find cells with highly confident label assignments through consensus and spreading these confident labels with a semi-supervised GCN. Using simulated data and two scRNA-seq datasets from different tissues, we show that our method improves accuracy over a simple consensus algorithm and the average of the underlying tools. We also compare our method to a nonparametric neighbor majority approach, showing comparable results. We then demonstrate that our GCN method allows for feature interpretation, identifying important genes for cell type classification. We present our completed pipeline, written in PyTorch, as an end-to-end tool for automating and interpreting the classification of scRNA-seq data. AVAILABILITY AND IMPLEMENTATION: Our code for conducting the experiments in this paper and using our model is available at https://github.com/lewinsohndp/scSHARP.

The origins and functional effects of postzygotic mutations throughout the human life span.

Postzygotic mutations (PZMs) begin to accrue in the human genome immediately after fertilization, but how and when PZMs affect development and lifetime health remain unclear. To study the origins and functional consequences of PZMs, we generated a multitissue atlas of PZMs spanning 54 tissue and cell types from 948 donors. Nearly half the variation in mutation burden among tissue samples can be explained by measured technical and biological effects, and 9% can be attributed to donor-specific effects. Through phylogenetic reconstruction of PZMs, we found that their type and predicted functional impact vary during prenatal development, across tissues, and through the germ cell life cycle. Thus, methods for interpreting effects across the body and the life span are needed to fully understand the consequences of genetic variants.

TAD Evolutionary and functional characterization reveals diversity in mammalian TAD boundary properties and function.

Topological associating domains (TADs) are self-interacting genomic units crucial for shaping gene regulation patterns. Despite their importance, the extent of their evolutionary conservation and its functional implications remain largely unknown. In this study, we generate Hi-C and ChIP-seq data and compare TAD organization across four primate and four rodent species, and characterize the genetic and epigenetic properties of TAD boundaries in correspondence to their evolutionary conservation. We find that only 14% of all human TAD boundaries are shared among all eight species (ultraconserved), while 15% are human-specific. Ultraconserved TAD boundaries have stronger insulation strength, CTCF binding, and enrichment of older retrotransposons, compared to species-specific boundaries. CRISPR-Cas9 knockouts of two ultraconserved boundaries in mouse models leads to tissue-specific gene expression changes and morphological phenotypes. Deletion of a human-specific boundary near the autism-related AUTS2 gene results in upregulation of this gene in neurons. Overall, our study provides pertinent TAD boundary evolutionary conservation annotations, and showcase the functional importance of TAD evolution.

Diverse monogenic subforms of human spermatogenic failure.

Non-obstructive azoospermia (NOA) is the most severe form of male infertility and typically incurable. Defining the genetic basis of NOA has proven challenging, and the most advanced classification of NOA subforms is not based on genetics, but simple description of testis histology. In this study, we exome-sequenced over 1000 clinically diagnosed NOA cases and identified a plausible recessive Mendelian cause in 20%. We find further support for 21 genes in a 2-stage burden test with 2072 cases and 11,587 fertile controls. The disrupted genes are primarily on the autosomes, enriched for undescribed human "knockouts", and, for the most part, have yet to be linked to a Mendelian trait. Integration with single-cell RNA sequencing data shows that azoospermia genes can be grouped into molecular subforms with synchronized expression patterns, and analogs of these subforms exist in mice. This analysis framework identifies groups of genes with known roles in spermatogenesis but also reveals unrecognized subforms, such as a set of genes expressed across mitotic divisions of differentiating spermatogonia. Our findings highlight NOA as an understudied Mendelian disorder and provide a conceptual structure for organizing the complex genetics of male infertility, which may provide a rational basis for disease classification.

SATINN: an automated neural network-based classification of testicular sections allows for high-throughput histopathology of mouse mutants.

MOTIVATION: The mammalian testis is a complex organ with a cellular composition that changes smoothly and cyclically in normal adults. While testis histology is already an invaluable tool for identifying and describing developmental differences in evolution and disease, methods for standardized, digital image analysis of testis are needed to expand the utility of this approach. RESULTS: We developed SATINN (Software for Analysis of Testis Images with Neural Networks), a multi-level framework for automated analysis of multiplexed immunofluorescence images from mouse testis. This approach uses residual learning to train convolutional neural networks (CNNs) to classify nuclei from seminiferous tubules into seven distinct cell types with an accuracy of 81.7%. These cell classifications are then used in a second-level tubule CNN, which places seminiferous tubules into one of 12 distinct tubule stages with 57.3% direct accuracy and 94.9% within ±1 stage. We further describe numerous cell- and tubule-level statistics that can be derived from wild-type testis. Finally, we demonstrate how the classifiers and derived statistics can be used to rapidly and precisely describe pathology by applying our methods to image data from two mutant mouse lines. Our results demonstrate the feasibility and potential of using computer-assisted analysis for testis histology, an area poised to evolve rapidly on the back of emerging, spatially resolved genomic and proteomic technologies. AVAILABILITY AND IMPLEMENTATION: The source code to reproduce the results described here and a SATINN standalone application with graphic-user interface are available from http://github.com/conradlab/SATINN. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Comparative single-cell analysis of biopsies clarifies pathogenic mechanisms in Klinefelter syndrome.

Klinefelter syndrome (KS), also known as 47, XXY, is characterized by a distinct set of physiological abnormalities, commonly including infertility. The molecular basis for Klinefelter-related infertility is still unclear, largely because of the cellular complexity of the testis and the intricate endocrine and paracrine signaling that regulates spermatogenesis. Here, we demonstrate an analysis framework for dissecting human testis pathology that uses comparative analysis of single-cell RNA-sequencing data from the biopsies of 12 human donors. By comparing donors from a range of ages and forms of infertility, we generate gene expression signatures that characterize normal testicular function and distinguish clinically distinct forms of male infertility. Unexpectedly, we identified a subpopulation of Sertoli cells within multiple individuals with KS that lack transcription from the XIST locus, and the consequence of this is increased X-linked gene expression compared to all other KS cell populations. By systematic assessment of known cell signaling pathways, we identify 72 pathways potentially active in testis, dozens of which appear upregulated in KS. Altogether our data support a model of pathogenic changes in interstitial cells cascading from loss of X inactivation in pubertal Sertoli cells and nominate dosage-sensitive factors secreted by Sertoli cells that may contribute to the process. Our findings demonstrate the value of comparative patient analysis in mapping genetic mechanisms of disease and identify an epigenetic phenomenon in KS Sertoli cells that may prove important for understanding causes of infertility and sex chromosome evolution.

Somatic PIK3R1 variation as a cause of vascular malformations and overgrowth.

PURPOSE: Somatic activating variants in the PI3K-AKT pathway cause vascular malformations with and without overgrowth. We previously reported an individual with capillary and lymphatic malformation harboring a pathogenic somatic variant in PIK3R1, which encodes three PI3K complex regulatory subunits. Here, we investigate PIK3R1 in a large cohort with vascular anomalies and identify an additional 16 individuals with somatic mosaic variants in PIK3R1. METHODS: Affected tissue from individuals with vascular lesions and overgrowth recruited from a multisite collaborative network was studied. Next-generation sequencing targeting coding regions of cell-signaling and cancer-associated genes was performed followed by assessment of variant pathogenicity. RESULTS: The phenotypic and variant spectrum associated with somatic variation in PIK3R1 is reported herein. Variants occurred in the inter-SH2 or N-terminal SH2 domains of all three PIK3R1 protein products. Phenotypic features overlapped those of the PIK3CA-related overgrowth spectrum (PROS). These overlapping features included mixed vascular malformations, sandal toe gap deformity with macrodactyly, lymphatic malformations, venous ectasias, and overgrowth of soft tissue or bone. CONCLUSION: Somatic PIK3R1 variants sharing attributes with cancer-associated variants cause complex vascular malformations and overgrowth. The PIK3R1-associated phenotypic spectrum overlaps with PROS. These data extend understanding of the diverse phenotypic spectrum attributable to genetic variation in the PI3K-AKT pathway.

Rare mutations in the complement regulatory gene CSMD1 are associated with male and female infertility.

Infertility in men and women is a complex genetic trait with shared biological bases between the sexes. Here, we perform a series of rare variant analyses across 73,185 women and men to identify genes that contribute to primary gonadal dysfunction. We report CSMD1, a complement regulatory protein on chromosome 8p23, as a strong candidate locus in both sexes. We show that CSMD1 is enriched at the germ-cell/somatic-cell interface in both male and female gonads. Csmd1-knockout males show increased rates of infertility with significantly increased complement C3 protein deposition in the testes, accompanied by severe histological degeneration. Knockout females show significant reduction in ovarian quality and breeding success, as well as mammary branching impairment. Double knockout of Csmd1 and C3 causes non-additive reduction in breeding success, suggesting that CSMD1 and the complement pathway play an important role in the normal postnatal development of the gonads in both sexes.

The Emerging Role of Chromatin Remodeling Factors in Female Pubertal Development.

To attain sexual competence, all mammalian species go through puberty, a maturational period during which body growth and development of secondary sexual characteristics occur. Puberty begins when the diurnal pulsatile gonadotropin-releasing hormone (GnRH) release from the hypothalamus increases for a prolonged period of time, driving the adenohypophysis to increase the pulsatile release of luteinizing hormone with diurnal periodicity. Increased pubertal GnRH secretion does not appear to be driven by inherent changes in GnRH neuronal activity; rather, it is induced by changes in transsynaptic and glial inputs to GnRH neurons. We now know that these changes involve a reduction in inhibitory transsynaptic inputs combined with increased transsynaptic and glial excitatory inputs to the GnRH neuronal network. Although the pubertal process is known to have a strong genetic component, during the last several years, epigenetics has been implicated as a significant regulatory mechanism through which GnRH release is first repressed before puberty and is involved later on during the increase in GnRH secretion that brings about the pubertal process. According to this concept, a central target of epigenetic regulation is the transcriptional machinery of neurons implicated in stimulating GnRH release. Here, we will briefly review the hormonal changes associated with the advent of female puberty and the role that excitatory transsynaptic inputs have in this process. In addition, we will examine the 3 major groups of epigenetic modifying enzymes expressed in the neuroendocrine hypothalamus, which was recently shown to be involved in pubertal development and progression.

Beyond sequence variation: assessment of copy number variation in adult glioblastoma through targeted tumor somatic profiling.

Glioblastoma is the most common primary malignancy of the adult central nervous system. Gliomagenesis involves a complex range of alterations, including sequence changes, copy number variations (CNVs), and epigenetic modifications, that have clinical implications for disease classification and prognosis. Thus, multiple testing modalities are required to support a complete diagnostic workup. The goal of this study was to streamline the multipart workflow by predicting both sequence changes and CNVs (specifically EGFR amplifications) from a single next-generation sequencing (NGS) test. Eighty-six primary and secondary glioblastomas were submitted for clinical NGS to report sequence variants from a concise panel of cancer-relevant genes. Most specimens underwent concomitant testing by methylation-specific polymerase chain reaction, immunohistochemistry, and fluorescence in situ hybridization. Using data generated during the course of clinical testing, we found that NGS-based variant predictions were concordant with immunohistochemistry and fluorescence in situ hybridization for IDH mutation and EGFR amplification status, respectively. We also noted that EGFR amplifications correlated with polysomy of chromosome 7, 19, and 20, and loss of PTEN and CDKN2A. EGFR-unamplified cases had lower rates of chromosome 7 polysomy, and PTEN and CDKN2A loss, but more CNVs overall. TP53, NF1, ATRX, and PDGFRA mutations were nearly exclusive to specimens without EGFR amplification. EGFR amplification was not associated with longer progression-free survival in this cohort, but amplifications were enriched in a group with slightly longer overall survival despite radiographic evidence of disease progression. Further study is needed to explore the mechanisms responsible for noted patterns of co-occurring variants and to correlate them with specific clinical outcomes.

Analyzing the Genetic Spectrum of Vascular Anomalies with Overgrowth via Cancer Genomics.

Vascular anomalies are variably associated with overgrowth, skeletal anomalies, and abnormalities of the brain, leptomeninges, and eye. We assembled a 16-institution network to determine the range of genetic variants associated with a spectrum of vascular anomalies with overgrowth, ranging from mild to severe. Because of the overlap between cancer-associated variants and previously described somatic variants in vascular overgrowth syndromes, we employed tumor genetic profiling via high-depth next-generation sequencing using a panel to assay affected tissue from a diverse cohort of subjects with vascular anomalies with overgrowth. Seventy-five percent (43/57) harbored pathogenic or likely pathogenic variants in 10 genes. We identified two genes (mTOR, PIK3R1) and several variants previously described in the setting of cancer but that, to our knowledge, have not been described in vascular malformations. All were identified at low variant allele frequency consistent with somatic mosaic etiology. By leveraging somatic variant detection technology typically applied to cancer in a cohort inclusive of broad phenotypic severity, we demonstrated that most vascular anomalies with overgrowth harbor postzygotic gain-of-function mutations in oncogenes. Furthermore, continued interrogation of oncogenes in benign developmental disorders could provide insight into fundamental mechanisms regulating cell growth.

FGFR2 amplification in colorectal adenocarcinoma.

FGFR2 is recurrently amplified in 5% of gastric cancers and 1%-4% of breast cancers; however, this molecular alteration has never been reported in a primary colorectal cancer specimen. Preclinical studies indicate that several FGFR tyrosine-kinase inhibitors (TKIs), such as AZD4547, have in vitro activity against the FGFR2-amplified colorectal cell line, NCI-H716. The efficacy of these inhibitors is currently under investigation in clinical trials for breast and gastric cancer. Thus, better characterizing colorectal tumors for FGFR2 amplification could identify a subset of patients who may benefit from FGFR TKI therapies. Here, we describe a novel FGFR2 amplification identified by clinical next-generation sequencing in a primary colorectal cancer. Further characterization of the tumor by immunohistochemistry showed neuroendocrine differentiation, similar to the reported properties of the NCI-H716 cell line. These findings demonstrate that the spectrum of potentially clinically actionable mutations detected by targeted clinical sequencing panels is not limited to only single-nucleotide polymorphisms and insertions/deletions but also to copy-number alterations.

Routine use of clinical exome-based next-generation sequencing for evaluation of patients with thrombotic microangiopathies.

Next-generation sequencing is increasingly used for clinical evaluation of patients presenting with thrombotic microangiopathies because it allows for simultaneous interrogation of multiple complement and coagulation pathway genes known to be associated with disease. However, the diagnostic yield is undefined in routine clinical practice. Historic studies relied on case-control cohorts, did not apply current guidelines for variant pathogenicity assessment, and used targeted gene enrichment combined with next-generation sequencing. A clinically enhanced exome, targeting ~54 Mb, was sequenced for 73 patients. Variant analysis and interpretation were performed on genes with biological relevance in thrombotic microangiopathy (C3,CD46, CFB, CFH, CFI, DGKE, and THBD). CFHR3-CFHR1 deletion status was also assessed using multiplex ligation-dependent probe amplification. Variants were classified using American College of Medical Genetics and Genomics guidelines. We identified 5 unique novel and 14 unique rare variants in 25% (18/73) of patients, including a total of 5 pathogenic, 4 likely pathogenic, and 15 variants of uncertain clinical significance. Nine patients had homozygous deletions in CFHR3-CFHR1. The diagnostic yield, defined as the presence of a pathogenic variant, likely pathogenic variant or homozygous deletion of CFHR3-CFHR1, was 25% for all patients tested. Variants of uncertain clinical significance were identified in 21% (15/73) of patients.These results illustrate the expected diagnositic yield in the setting of thrombotic microangiopathies through the application of standardized variant interpretation, and highlight the utility of such an approach. Sequencing a clinically enhanced exome to enable targeted, disease-specific variant analysis is a viable approach. The moderate rate of variants of uncertain clinical significance highlights the paucity of data surrounding the variants in our cohort and illustrates the need for expanded variant curation resources to aid in thrombotic microangiopathy-related disease variant classification.

Targeted Next-Generation Sequencing in Molecular Subtyping of Lower-Grade Diffuse Gliomas: Application of the World Health Organization's 2016 Revised Criteria for Central Nervous System Tumors.

The 2007 World Health Organization Classification of Tumours of the Central Nervous System classifies lower-grade gliomas [LGGs (grades II to III diffuse gliomas)] morphologically as astrocytomas or oligodendrogliomas, and tumors with unclear ambiguous morphology as oligoastrocytomas. The World Health Organization's newly released (2016) classification incorporates molecular data. A single, targeted next-generation sequencing (NGS) panel was used for detecting single-nucleotide variation and copy number variation in 50 LGG cases originally classified using the 2007 criteria, including 36 oligoastrocytomas, 11 oligodendrogliomas, 2 astrocytomas, and 1 LGG not otherwise specified. NGS results were compared with those from IHC analysis and fluorescence in situ hybridization to assess concordance and to categorize the tumors according to the 2016 criteria. NGS results were concordant with those from IHC analysis in all cases. In 3 cases, NGS was superior to fluorescence in situ hybridization in distinguishing segmental chromosomal losses from whole-arm deletions. The NGS approach was effective in reclassifying 36 oligoastrocytomas as 30 astrocytomas (20 IDH1/2 mutant and 10 IDH1/2 wild type) and 6 oligodendrogliomas, and 1 oligodendroglioma as an astrocytoma (IDH1/2 mutant). Here we show that a single, targeted NGS assay can serve as the sole testing modality for categorizing LGG according to the World Health Organization's 2016 diagnostic scheme. This modality affords greater accuracy and efficiency while reducing specimen tissue requirements compared with multimodal approaches.

A protein allergen microarray detects specific IgE to pollen surface, cytoplasmic, and commercial allergen extracts.

BACKGROUND: Current diagnostics for allergies, such as skin prick and radioallergosorbent tests, do not allow for inexpensive, high-throughput screening of patients. Additionally, extracts used in these methods are made from washed pollen that lacks pollen surface materials that may contain allergens. METHODOLOGY/PRINCIPAL FINDINGS: We sought to develop a high-throughput assay to rapidly measure allergen-specific IgE in sera and to explore the relative allergenicity of different pollen fractions (i.e. surface, cytoplasmic, commercial extracts). To do this, we generated a protein microarray containing surface, cytoplasmic, and commercial extracts from 22 pollen species, commercial extracts from nine non-pollen allergens, and five recombinant allergenic proteins. Pollen surface and cytoplasmic fractions were prepared by extraction into organic solvents and aqueous buffers, respectively. Arrays were incubated with <25 uL of serum from 176 individuals and bound IgE was detected by indirect immunofluorescence, providing a high-throughput measurement of IgE. We demonstrated that the allergen microarray is a reproducible method to measure allergen-specific IgE in small amounts of sera. Using this tool, we demonstrated that specific IgE clusters according to the phylogeny of the allergen source. We also showed that the pollen surface, which has been largely overlooked in the past, contained potent allergens. Although, as a class, cytoplasmic fractions obtained by our pulverization/precipitation method were comparable to commercial extracts, many individual allergens showed significant differences. CONCLUSIONS/SIGNIFICANCE: These results support the hypothesis that protein microarray technology is a useful tool for both research and in the clinic. It could provide a more efficient and less painful alternative to traditionally used skin prick tests, making it economically feasible to compare allergen sensitivity of different populations, monitor individual responses over time, and facilitate genetic studies on pollen allergy.